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Abstract. Wc present a new path integral method to analyze stochastically perturbed ordinary 
differential equations with multiple time scales. The objective of this method is to derive from 
the original system a new stochastic differential equation describing the system's evolution on slow 
time scales. For this purpose, wc start from the corresponding path integral representation of the 
stochastic system and apply a multi-scale expansion to the associated path integral kernel of the 
corresponding Lagrangian. As a concrete example, we apply tins expansion to a system that arises 
in the study of random dispersion fluctuations in dispersion-managed fiber optic communications. 
Moreover, we show that, for this particular example, the new path integration method yields the 
same result at leading order as an asymptotic expansion of the associated Fokker-Planck equation. 



1. Introduction. Physical phenomena exhibiting scale separation are ubiqui- 
tous in a wide range of fields, most obviously in those permitting a description via 
continuum mechanics such as fluid dynamics, electromagnetism, and material science. 
The ability to separate the (macroscopic) evolution on slow time scales from the (mi- 
croscopic) evolution on fast time scales is key to understanding a system's behavior 
or to performing numerical simulations in an efficient way. Many mathematical tech- 
niques have been developed for explicitly performing this separation of scales in the 
absence of stochastic fluctuations, e.g. multi-scale expansions [1], Lie transform [2], 
and renormalization group methods [3]. In a typical multi-scale analysis, we start 
from a given system, which we write here as an ordinary differential equation for a 
vector X = x{t), 

where e is a small parameter. Upon introducing, for instance, two time scales to = t 
and ti = e^t (the actual power of 7 > is selected by the problem), we derive 
systematically a new ordinary differential equation 

describing the evolution of the system on the slow time scale represented by ti through 
X = X{ti). In this article, we are interested in the effect of stochastic perturbations 
on this description over the coarse-grained (slow) scales. Intuitively we expect to 
derive from the microscopic SDE given by 

dx 

— = f{x, t, e) + g{x, t, e)^{t) , 
where ^ represents white noise, a new stochastic differential equation of the form 

^^=F{X,h) + G{X,h)E{h), 
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describing the slow scale evolution, where S represents white noise on the slow time 
scale ti. Our intent is not to explore under what circumstances or over what range 
of parameters such an approximation is valid; rather, we present a practical method 
based on path integrals that can be used for explicitly calculating the coarse-grained 
stochastic differential equation from the original system, implicitly assuming that such 
an approach is valid. 

For this purpose, we now state the problem in a slightly more formal way: Assume 
a system of stochastic ordinary differential (Langevin) equations expressed by 

f{x,t;e) + V~eg{x,t)m, x{0) ^ a (1.1) 

where x € M", t G M+, / : [M" x M+] W\ g : [M" x M+] and / is 

assumed to be sufficiently differentiable in e. The random term, ^(i), is assumed to 
be m-dimensional delta-correlated white noise, i.e., with 

< UW) >^ Sit ~ (1-2) 

where < • > denotes the ensemble average, and Gaussian-distributed with a proba- 
bility density given by 

PiO-j^^.eM'l^-O- (1.3) 

We choose to take the Stratonovich interpretation of Eqn. |1.1| given our interest in 
modeling physical processes; however, this is simply a convention given the general 
form of Eqn. |1.1[ Further assumptions will be indicated as they become necessary. 

The goal is to obtain another stochastic ODE that captures the evolution of x on 
a slow scale when < e <JC 1, i.e., to "scale up" the noise from its microscopic rep- 
resentation to a macroscopic noise that reflects the fast microscopic response of x in 
addition to the noise itself. The ability to integrate a macroscopic stochastic equation 
on a much slower time scale, where the microscopic trajectories can be reconstructed 
approximately from the macroscopic trajectories, is of obvious benefit from the per- 
spectives of asymptotic analysis and numerical generation of the statistical properties 
of X. 

For the reader who is less acquainted with the scaling of small parameters in 



stochastic equations, we comment briefly on the factor ^/e in the problem ( 1.1 1 in the 
context of an expansion using multiple time scales. A naive expansion of x in Eqn.fTT] 
in powers of e (a/c) is quickly seen to become disordered when t ^ 1/e {t ^ l/v^)- 
A standard method that is widely used in the deterministic case [1] to address this 
onset of secularity is the use of multiple scales in time, i.e., let 

{to,h,...) + ex^^\to,ti,...) + (1.4) 



X 



(0) 



where tj — eH. Whereas the dependence of x^-^\ x^-^\ etc., on the "fast" time variable 
to alone is insufficient to prevent this expansion of x from becoming disordered, their 
dependence on the "slow" time variables ti, etc., allows sufficient freedom to suppress 
secular growth up to times of order 1/e. Before proceeding, we have to justify the 



particular expansion suggested in Eqn. 1.4 For simplicity, we do this in a one- 



dimensional context; the argument generalizes in an obvious way to higher dimensions. 
Were ^{t) a deterministic forcing term, we would ordinarily expand x and t in powers 
of ^/e. In this case, however, ^(i) is formally the time-derivative of a Wiener process 
W{t), with Eqn. [lT] more appropriately expressed in the form 

dx = f{x,t;e)dt + ^/lg{x,t)dW, (1.5) 
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where {W{t)W{t')) = min(t,t'), i-e., with {dW"^) = {{W{dt) - W{{))Y) = dt. As x{t) 
is a stochastic process, we wish to track all of its moments, as an alternative to the 
direct evolution of the probability density function via the Fokker-Planck equation 
(see Sec. |4|. In a time step of dt, we have 

{x{t + dt) - x{t)) = (dx) ^ {f{x, t; e)) dt, and 
{x'^{t + dt) -x^{t)) = 2{xdx) + {dx^) 

^2{xf{x,t;e))dt + e{g^{x,t))dt, 

so that all terms of 0{e^^'^) vanish for j odd. An expansion in powers of e is therefore 
self-consistent . We also see from this consideration that the influence of the stochastic 
driving term in Eqn. |l.l| is first felt at 0(e), suggesting that the term of 0{^/edW) 
mimics a term of 0(edt). 

The paper is organized as follows: In Sec. [2] we present the path integral method 
for deriving the slow stochastic dynamics from the original system. In Sec. [3] we 
apply this method to a particular problem in fiber optics. In what follows, we show 
that the new method, for the chosen example, yields the same result as an asymp- 
totic expansion of the associated Fokker-Planck equation (Sec. Conclusions are 
presented in Sec. [5] 

2. Averaging of path integral kernels. Instead of considering the stochastic 
differential Eqn. |1.1| directly we can consider the probability density p{x, t) corre- 
sponding to the stochastic process x and write down a path integral representation 
for p{x,t) [4-6]: 

p{x,t)^ [ l?2;(r) e--/'o = /" Vx{t)V{x{t)) (2.1) 

JCix,t\a,0) Jc{x,t\a,0) 

with the Lagrangian L given by 

^ = ^ E ^"'^ {-^ {-f^ + "Sf^'^) ) ' ' ('-'^ 

written according to Einstein's summation notation (see Appendix). Here we set s = 
1/2 according to the Stratonovich interpretation of Eqn. |l.l[ In cases where h = 
does not exist, it is necessary to consider a regularized path integral representation. 
A short derivation for both cases can be found in the appendix. In Eqn. |2.1[ we take 
all paths r e C{x,t\a,0) into account, that lead from the initial point (a, 0) to the 
final point {x,t). We want to develop a multi-scale technique on the level of this 
path integral representation and assume a periodicity on the fast time scale to with 
a period t* . Assume that the final time t — {K -\- Vjt* , meaning that we propagate 
for K + 1 periods and with K intermediate times tofc = kt* . The fast time scale 
to = t characterizes the system's evolution within one period, e.g. for a time interval 
[tok,tok + t*]. The evolution on the time scale ti = et characterizes the slow evolution 
of the system. For a given path F g C{x, t\a, 0), we take Xk — Xk = x{kt*), < k < K 
as sample points of the process x at times ti^ = ekt* and we are looking for a slowly 
varying process X = X(ti) as a continuum limit of this sampled process (Xk). We 
can write down the transition probability p{x,t) = p(x,t\a,0) as a X-dimensional 
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integral 



ith 



p{x,t\a,0) = J + l)t*\xk,kt*)^ dxi...dxK 

= J (^Ylp{Xk+i,tnk+i)\Xk,tik)^ dXi...dXK 

p{Xk+i,ti(k+i)\Xk,tik) = p{xk+i, {k + l)t*\xk,kt*). 



(2.3) 



In order to use the path integral representation given in Eqn. 2.1 we can write any 
path r G C{x, t\a, 0) as 



K 



r = ^r,, Tkec{xk+i,{k + i)t*\xk,kt*). 



k=0 



The following Figure 2.1 illustrates this decomposition of a path into polygons con- 
necting the start- and endpoints of the spans and the paths within one period. The 




kf (k+l)t* t 

Fig. 2.1. A path from (a, 0) to {x,t) can be decomposed into a polygon connecting the 
beginning and the endpoints of the periods and paths within the periods. The short-ttme 
propagator averages over all paths within one period. As an example, the paths from 
time kt* to time (k + l)t* are drawn as solid lines. 

transition probability p^x^^i, (fc + l)t* \xk, kt*) from the beginning of a period to its 
end is then given by the path integral 

p{xk+i,ik+l)t*\xk,kt*)^ f V{x{t))Vx{t) (2.4) 

Jc(xk + i,(k+l)t*\xkM') 
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and in the continuum limit for the slow scale ti, the polygons connecting (a, 0) with 
{x,t) through K intermediate points become paths X{ti) and the probabiHty distri- 
bution for the polygons becomes the probability distribution V{X{ti)) of these paths 



^ (f[piXk+i,t^k+i)\Xk,hk)]^P{X{h)). (2.5) 

K — »oo. — n) \ I 
\k=l I 

Thus we obtain the final representation of the transition probability on the slow scales: 
p(x,t)=p(X,ii|a,0)= / V{X{Tr))VX{Tr) (2.6) 

JC{X,t^\afi) 

The last three equations are the mathematical formulation of the intuitive approach: 
First, one considers the dynamics of the system on the fast scales. The resulting path 
integral yields the propagator of the system from the beginning of a period to the end 
of the period. Combining these propagators leads to the probability of the path on 
the slow scales. Therefore, as an abbreviation for combining Eqn. |2.4| and Eqn. |2.6[ 
we can write the transition probability p as a hierarchy of path integrals: 



p(a;,i|a, 0) 



/ (l™n / V{X(t^))VX{t^)\^ VX{n) (2.7) 



This method can be generalized in an obvious way: Introducing J time scales tj = e^t, 
the product of the propagators on the time scale tj leads to the probability of a path 



on the time scale ij+i. Using the same abbreviated notation as in Eqn. 2.7 
represent p = p{x, t\a, 0) as 



P = 



J ... (^limj] I (liniH / ViX{To))VXiTo)^ 2?^(ti)) -VXiT,) (2.8) 



Once we have obtained the path integral representation on the slow scale with an 
averaged kernel 

P{X{ti)) = e^-fo' i(X(ri),X(ri),ri)dri (2.9) 

we can also convert back in order to find the corresponding stochastic equation on 
the slow time scale. 

3. Dispersion-managed fiber-optic communications. To demonstrate the 
application of the method described above, we consider the propagation of an electro- 
magnetic field through an optical fiber with a periodic, piecewise-constant coefficient 
of dispersion. The model equation for this process is the dispersion-managed nonlinear 
Schrodinger (DMNLS) equation given by 

iE, + d{z)Ett + e\E^\E^O, (3.1) 

where z, t, E{t,z), d{z) are dimensionless quantities representing distance down the 
fiber (the "time-like" evolution variable in signaling coordinates), time (the trans- 
verse variable in signaling coordinates), the electric field envelope and the dispersion 
coefficient. The scaling parameter < e <C 1 reflects the typical case of a weak 
nonlinearity relative to the dispersion strength at any point in the flber. Under typ- 
ical deterministic dispersion management, the fiber dispersion has a small positive 
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(anomalous) mean value and 0(1) local mean-zero value, i.e., d{z) — ed^v + do{z), 
where, for instance. 



-d 0<z<l/4 
do{z) ^ {-d 1/4 < z < 3/4 , (3.2) 
-d 3/4 < z < 1 

and where 1^0(2) has unit periodicity. Here, we consider the case where d{z) has an 
additional random component, giving d{z) ~ edav + do{z) + V eD^{z), where ^(z) is 
taken to be (5-correlated white noise with unit strength, i.e., 

<e(z)e(^')>-'5(z-z'). (3.3) 

For more details on the derivation of this model from Maxwell's equations, refer 
to [7,8]. It has already been demonstrated [9,10] that the effect of this random 
dispersion on single periodic solutions (referred to as "dispersion-managed solitons" ) 
of Eqn. |3.1| can effectively be captured by the two-dimensional system of stochastic 
ordinary differential equations given by 

dT 

— = 4d(z)M and 
dz 

dM _ Cid(z) eCa 
~ 4r3 4r2 ' 

where T(z) and M(z) represent the width and chirp (i.e., the quadratic component of 
the phase) of the soliton. For derivations of this system using a variational principle, 
by applying the lens transform followed by expansion in Gauss-Hermite modes, or 
using moments with a closure condition, see [11, 12]. Constants Ci and C2 depend 
on the particular shape assumed for the soliton. When e = 0, these equations sim- 
ply reflect the linear Schrodinger equation with mean-zero dispersion, and all initial 
conditions, in particular those of the form {r(0), M(0)} = {Tq, 0}, yield periodic so- 
lutions. When e is finite and D = (i.e., the dispersion is deterministic and periodic), 
the nonlinearity and finite-mean dispersion "select" a member of this one-dimensional 
family of periodic solutions that persists for e > [13]. When e > and D — 0(1), 
the width and chirp of this soliton undergo a random walk in addition to the "breath- 
ing" induced by the deterministic component of the dispersion map. We now seek to 
apply the methods discussed in the previous sections to determine the properties of 
this random walk. 

Before we proceed, it is useful to note that the action-angle coordinates [14] in 
the e = case also simplify analysis for the finite-e case. Letting = Ci/T'^ + 16M^ 
and [3 = ATM gives 

dfl ^ dh 

~T~ = and 
dz op 

dp , , s,^ dh 

£ = dom-^g^ (3.4) 



where 



2C n^/"^ 1 

'^("'^)-3(crw^-2^-"'- ^'-'^ 
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An initial condition in Eqns. 3.4 equal to {TojO} corresponds to an initial condition 
in Eqns.jsljof {f^o,0}, where = Ci/Tg. 

Intuitively, we expect the solutions of the stochastic dynamical system presented 



by (3.4 1 to experience scale-separation as in the deterministic case. This separation 



of scales should result in fast variations of the probability densities for /3 and Q 
within a span and a slow evolution of both quantities over many periods visible in the 
corresponding Poincare sections. Numerical simulations confirm this intuition. Fig. 
3.1|and fig. |3.2| present the evolution of both probability densities. 



1.5 
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Fig. 3.1. Evolution of the probability density of fl 



We start our m ulti- scale analysis of the system by immediately writing Eqns. 3.4 
in the form of Eqn. with x = {ft (3)'^ and t = z. Thus, 

(3.7) 



^) ^ W{z)n ) + ' [-dh/dn ) ' ^^'^ 



g{x,z) 





n 



where we write £,{z) = {^i{z) £,2{z))'^ formally even though only ^2 is relevant. The 
initial condition is simply a — {^Iq 0)'^. Note that the behavior at leading order is 
easily solved to yield 



,(0) 



1 

Rizo) 



and 



i?(zo) 



doiO dC- 



(3.8) 



In order to obtain coarse-grained equations from the original system given by 
Eqn.|3.4"l we first express the transition probability in the form of a path integral, i.e.. 



p{x,z) = / V(3{0 e-A^o MP-MOn+.mTdc ^ ( _ _ /'s^l^dC 

Jc[f3,z\l3o,0] \ Jo 

(3.9) 



-5 
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Fig. 3.2. Evolution of the probability density of [3 



with the function h given by Eqn. |3.5| A short derivation of this representation is 
found in the appendix. For the short-time propagator we note in particular that, since 

( r''+^ dh \ 

X 5in{k+l)~n{k)~ J 3e— dCj. (3.10) 

We wish to apply a semi-classical approximation to this short-term propagator by 
finding the minimal trajectory using the associated Euler-Lagrange equations [5] . For 
this purpose, we consider a multi-scale expansion 

n = r!("'(^o, zi) + efl'^^Hzo, zi) + ... (3.11) 
(3 = /3(°)(zo, zi) + ep^'Hzo, z,) + . . . (3.12) 

and find at the leading order 

n^°'>{zo,zi)^n{zi) and P^'^'' {zq, zi) = R{zo)n{zi) + ^{zi) . (3.13) 

Therefore, within one period, can be considered to be constant and the constraint 
fl — 3edh/dp yields an evolution equation for the zi dependence 

n'.f=r3|.=, (3.14) 

dzi d(3 

as well as the solution of the first-order correction fi*-^^ 

f](i)(zo,zi)- / 3—dI-n'{z,)za. (3.15) 
Jk op 



Here, we have applied standard deterministic multi-scale analysis - consistent with 
the fact that the evolution equation of does not contain any further source of noise. 
We can now rewrite the Lagrangian in the kernel of the above path integral given in 
Eqn. 3.10| by noting that 



In this way, the Lagrangian depends on /^^^HC) and (if' = dp^^^/dC,, or more formally 
written as 



and the corresponding Euler-Lagrange equation is given by 

d dL dL 



Since (3^^^^ is cyclic, we immediately find 

5/?(i) /■'' dh - dh 

m' + p' + ^^ d{c) 3^d~z + dion'c + ^ = c (a.i?) 

The value of the constant C is found by the boundary conditions on /J'^) and for 
/3(i)(fc + 1, zi) = P^^\k, zi) we find by integration of the above equation over the fast 
variable s immediately 

^-^'+/Y3i + gVc (3.18) 



df3 an 

which finally lets the averaged Lagrangian become 



- _ 1 1 



-, , '=+1 f dh dh , ,^ 



1 



(3.19) 



This, however corresponds again exactly to a stochastic equation with a noise process 
^ — (^{z\) on the slow scale given by 

and together with Eqn. |3.15[ we have found the coarse-grained system of stochastic 
differential equations. Note that in the present case, it is possible to solve all integrals 
analytically such that we can rewrite the system in terms of the slow variable y with 
= f2 and 2/2 = /3 as 



A {v^\ ^ /^ri(yi>2/2)A ^ / 

9 



(3.21) 



The slow-scale white noise process S is simply given by 

The functions Fi and r2 can be expressed as 
6 



ri(2/) 



dyi 
d 



dH, 1 ^ , 



-Q^{yi,y.--^dy,) 



Ci + (2/2 + -Jyif 



-1/2 



(3.22) 



(3.23) 



(3.24) 



and 

My) 



3 

M 
3 

d^yi 



H{yi,y2 + ^dyi) - H{yi,y2 - ^dyi) 



^(j/1,2/2 + ^dyi) + —{yi,y2- -dyi) 



davyi 



(3.25) 



sinh 



-1 f y2 + \dyi 



sinh 



_i / 2/2 - jdy i 



-C2^fy\' 



Cx + (?/2 + ^rfyi)^ 



-1-1/2 



+ 



Cx + (y2 - \dyxf 



1-1/2- 



with 



i/(f^,/3) = -C2f^3/2sinh-i(/3/ycr), 



(3.26) 



(3.27) 



Although the above analysis was carried out for a particular example, it shows the 
steps that are necessary to apply and modify the method for other cases. Note that it 
is not a requirement that the occurring integrals can be solved analytically. Moreover, 
using the semi-classical method in order to calculate the short-term propagator is 
appropriate in this case since the evolution of the system is linear on a short time 
scale. For systems that show nonlinear behavior on short times, other methods of 
evaluating the corresponding path integral can be applied. 

4. Perturbative analysis of Fokker-Planck equation. In this section we 
show that a perturbative analysis of the corresponding Fokker-Planck equation will 
yield for the particular example exactly the same result as the path integral based 
method. The idea to perform asymptotic analysis on the level of the Fokker-Planck 
equation in order to characterize the behavior of the associated stochastic process has 
been applied very successfully in a variety of contexts [15]. Since the Fokker-Planck 
equation is deterministic, we can use the classic theory of deterministic multi-scale 
expansions, adapted to the Fokker-Planck equation. 

Hence, instead of considering the path integral representation of the stochastic 
process a;(t), we now are considering the evolution equation of the probability density 
•p(x, t) given by the Fokker-Planck equation 



dp d_ 1 d 

dt~ dxi^^'^' ^ 2^ dxi 
10 



(4.1) 



with initial condition p{x,0) — Y[i=i^i^i ~ Note that, as before, we are using 
Einstein's summation convention and the Stratonovich interpretation of Eqn. [O] [16], 
Again, the secular growth we obtain when trying a naive perturbation expansion 
for p suggests we adopt a multiple-scales approach: 



p{x,t)=p'^"\x,to,h,...)+ep'^'>{x,ta,ti,...) + 



(4.2) 



where tj — cH. In this case, we have no ambiguity in the scale at which the effect 
of randomness is felt, since the Fokker-Planck equation is completely deterministic. 
The probability density p captures all moments of x(t). The leading order of this 
expansion satisfies 



dp(0) 



d 



dtQ dx 
so that along characteristics x — Xc{to; y) where 
dx. 



Mx,to;0)p^"^ 



we have that 



dto 

dp^ 
dto 



f{xc{to),to;0) and Xc{0) = y, 



dx, 



MxcitQ),tQ;0) 



P 



(0) 



(4.3) 



(4.4) 



(4.5) 



The tg dependence has been written as a full derivative here to denote the fact that 
it is taken along characteristics; naturally, p also depends on the slower time scale ti. 
The solution is then 



=P(2/,ii,...)exp 

with initial condition 



_d_ 

dx, 



/,(x,(c),C;0) dc 



(4.6) 



P{y,0) = l[5{y,-a,). 



(4.7) 



The next order is 



Lp 



(1) - 



where 



with adjoint operator 



_d_ 
dxi 



L = 



it 



de 



€=0 



l_d_ 
2 dxj 



d_ d_ 
dto di 



[Mx,t;0)o] 



d d 
— - Mx,t;0) — 
dto dxi 



d 



(4.8) 



(4.9) 



(4.10) 



By comparison with the leading order, it is clear that ker is spanned by functions 
that are constant in tp along characteristics xdto; y), i.e., any square-integrable func- 
tion (t>{y) (in fact, we have thus far omitted a discussion of the appropriate boundary 
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conditions, but for simplicity wc will assume i*-pcriodicity in and sufficiently fast 
decay in |x|). The Fredholm alternative therefore requires that 



= 



dP 



exp 







dx. 



■/,(x,(C),C;0)dC 



d_ 
dxi 

l_d_ 



df, 
de 



_d_ 

dXn 



P{y,ti)exp 
'ik — i9jkPiy,ti)) exp 



/,(x,(C),C;0)dC 



' dxi 







dxi 



fiix,iOX;0)dC 



(4.11) 



for all members (j>{y) of the admissible class of functions. This clearly leads to an IB VP 
for P(y,ti) which is now, depending on the nonempty elements of the matrix g(x,t), 
parabolic rather than hyperbolic. The particular form of this equation is seen most 
easily when the characteristic coordinate transformation x = Xc(to; y) is invertible, as 
we can then write down directly the transformed equation in y and ti. This will be 
the case for our current example, as we will see below. 

In application the concrete problem given by Eqns. [3^ we note immediately from 
Eqn.lS^that 



d_ 

dxi 



Mx,z;0)^0, 



(4.12) 



which greatly simplifies matters. As before, the leading order characteristics are given 
by 



Xc{zo) = 



1 







y, 



(4.13) 



,^(^o) 1, 

with the leading-order solution simply given by 

p(o)(x,zo,zi,...)-P(y,zi,...). (4.14) 
After using the above characteristics to change variables from x to y, the dependence 



of P on y and zi, as imposed by the Fredholm alternative condition in Eqn. 4.11 
given by 



IS 



dP 

dzi 



dzo 
d 

dyi 



oyi aV2 



dh 
dp 



{yi,V2 + R{zo)yi)P{y-,zi) 



dh 



{yi,y2 + R{zo)yi)P{y,zi) 



1 ^(Pf>_ 



dP d 

-Ctavyi -7. 3 



oy2 oyi Jo op 



_d 

dy2 Jo 



(•,.) + 3i?(zo)^(-,-) 



P{y, zi) dzo 



1 f)'^P 



where the derivatives of H{Q.,(3) are evaluated along the characteristics, that can be 
written explicitly as O = j/i, /3 = y2 + R{zo)yi- Finally, after some algebra we observe 
that the above evaluates to 



dP 
dzi 



oyi dy2 2 dy^ 

12 



(4.16) 



which is exactly the Fokker-Planck equation for the slow process given by Eq. |3.21| 
obtained in Sec.[3l 

5. Conclusion. We have presented a new method to derive equations that de- 
scribe the slow evolution of a system of stochastic ordinary differential equations. The 
conceptual basis of the method is the expression of state transition probabilities as 
path integrals, where a separation of scales allows a decomposition of the paths into 
self-similar segments. By performing a multi-scale expansion of the Lagrangian in 
each segment, we obtain an effective path integral that represents the coarse-grained 
noise process relevant to the slow dynamics of the system. From the appropriate limit 
of this path integral (i.e., refinement of the path partition), we obtain a new system of 
stochastic differential equations that describe behavior over long time scales. We have 
applied our method to a stochastic system arising in nonlinear optics, and have shown 
that our method produces a leading-order result that is identical to that obtained from 
a standard asymptotic analysis of the singularly perturbed Fokker-Planck equation 
associated with the stochastic system. The conditions necessary for this agreement 
to hold is a topic of ongoing research. 
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Appendix A. Path integral representation. 

In this appendix we give a brief overview of how to derive a path integral repre- 
sentation of a stochastic system following [17]. We start with a stochastic equation of 



the form ( 1 . 1 1 , 



x = f{x,t;e) + ^g{x,t)^{t), x(0) = a, (A.l) 
where we have denoted the noise strength by rj instead of e to reflect the fact that 



noise need not be small in this formalism. We take (A.l I to be the limit of difference 



equations for Axl*^^ (where the subscript and superscript now represent the component 
index and the time increment, respectively) where 



Aa; 



(fe) 



dg.. 



(fc) \ 



dxi 



(A.2) 



The time increments At identify a regular partition of interval [0, t] with K inter- 
mediate points. All terms in the stochastic difference equation are evaluated at the 
pre-point, hence 

and the parameter s accounts for the different possible definitions of the stochastic 
integral: a choice of s = corresponds to the Ito interpretation of ( A.l ) and a choice 
of s = 1/2 corresponds to the Stratonovich interpretation. We can rewrite Ax) 

(k) 

using the fact that the Brownian increments Aw^ are uncorrelated and of order 
VAt. This yields 



A (fc) /^-f(fc) , (fc)\ A + 



r- (fc) A (fc) 



Defining now the inverse of g to be ft, = g ^ , we can solve for the Aw, 

(fc) / / aJfc) 



(fc) 



(A.3) 



and obtain 



Aw 



(fc) 




+ dxi 



At 



(A.4) 



Note that if g is not invertible, an appropriate regularization has to be chosen (see 
Appendix B). In this section, however, we will simply take g to be nonsingular for all 
times t and all points x. 

The joint probability distribution of a path [w) = {'w^^\ w'^^^-*) is given by 



K+l 



p(H) = n 



1 



(A.5) 



kA V2^f=i 

The derivation of a path integral representation for the stochastic differential equa- 



tion (A.l I amounts to performing a coordinate transformation in the above probability 
density from the path {w) to the path {x) given by [x) = {x'^^\ x^^+^^). The de- 



terminant I J| of the Jacobian of this transformation is simple, since the discretization 
is based on the pre-point. We find 



K+l 

1^1 =n 



dxi' 



dwi 



dx). 



K+l 



n ^-"/^det([ft(f ]). 



(A.6) 



fc=i 



14 



We now define the action S as the following limit for K ^ oo and At 0: 

K+l n 



S= lim y y (Au;f 

fc=l i=l 

n //,(ni) 

= lini y y 

k=\ i=l \ 

1 -f^+l " 9 / A 

E E ('.!?)" ^ 

' k=l i=l \ 



,(fe) 



= lim 



At 



2?7 

where we have suppressed the arguments {x{t)^t) in hik, fk, etc., for simplicity of 
presentation. Defining 



K+l 



Vx{t)^ lim TT ^-—^^dei{[h^if])f\dx. 



t ^ n 

n „ det(ffe„(x(T),T)])]Jda;,(T) 
r=0 f v27rdTj 



we can write the path integral representation for the transition probability from the 
initial condition .x(0) = a to the point x{t) ~ x as 



p{x,t) 



I Vx{t) e~ 

Jcix,t\a,0) 



L{x(r),x{T),T) 



(A.7) 



where C(x, t\a, 0) denotes the set of all possible paths that connect the point (a, 0) to 
the point (x, t). 

Appendix B. Regularized path integral representation 

of the (r2,/3) equations. In this section we show how to derive a path integral 
representation in (ri,/3) for the particular system given by Eqns. (3.4 1. For simplicity 
we rewrite these equations as 



dp 
dz 



(B.l) 
(B.2) 



where we obtain (3.4) for k = 0. Similar regularization techniques are well-established, 
e.g., in the case of Brownian motion [5]. The matrices g and h in (A. 4 1 are given by 

,-1 



9 



K 

n 



K - 

n-^ 



such that 



hik^ — gij = 

OXi 
15 



(B.3) 



(B.4) 



for i — 1,2. With x = {n,/3)'^ and a — (rioj/3o)"^, we can write the corresponding 
path integral representation of (B.l[ B.2| as 



c[x,z\a,o] H \/27rdC V27rdC 



In order to take the hmit k — > 0, we note that 



hm I I , _ 



(B.5) 



(B.6) 



as can been seen by straightforward calculations using the same discretization as in 
the derivation of the path integral representation. Integration over dfl{() yields then 



or, written in a more compact way. 



p{x,z) 



C[fj,z\f3o,0] 



(B.7) 
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